Quantum heating of a parametrically modulated oscillator: spectral signatures 
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We show that the noise spectrum of a parametrically excited nonlinear oscillator can display a fine 
structure. It emerges from the interplay of the nonequidistance of the oscillator quasienergy levels 
and quantum heating that accompanies relaxation. The heating leads to a finite-width distribution 
over the quasienergy, or Floquet states even for zero temperature of the thermal reservoir coupled 
to the oscillator. The fine structure is due to transitions from different quasienergy levels, and 
thus it provides a sensitive tool for studying the distribution. For larger damping, where the fine 
structure is smeared out, quantum heating can be detected from the characteristic double-peak 
structure of the spectrum, which results from transitions accompanied by the increase or decrease 
of the quasienergy. 



PACS numbers: 03.65.Yz,05.40.-a,42.65.-k,85.25.Cp 

I. INTRODUCTION 

Nonlinearity is advantageous for observing quantum 
effects in vibrational systems. It makes the energy levels 
nonequidistant and the frequencies of different inter-level 
transitions different, which in turn enables spectroscopic 
observation of the quantum energy levels. In addition, 
nonlinearity leads to an interesting behavior of vibra- 
tional systems in external periodic fields, including the 
onset of bistability of forced vibrations. The interest in 
quantum effects in modulated nonlinear oscillators sig- 
nificantly increased recently due to the development of 
high-quality microwave resonators with the anharmonic- 
ity provided by Josephson junctions and to applications 
of these systems in quantum information [THB]. The 
long-sought [7] quantum regime has been reached also in 
nanomechanical systems |S]-[TUj. This development has 
opened the possibility of measurements of a single quan- 
tum nonlinear oscillator rather than of an ensemble of 
oscillators. In addition, the systems are versatile and 
allow accessing different dynamical regimes. 

One of the important problems that can be addressed 
with strongly modulated nonlinear oscillators is quantum 
fluctuations far from thermal equilibrium. In addition 
to the standard quantum uncertainty, such fluctuations 
come from the coupling of a quantum system to a thermal 
bath. The coupling leads to relaxation of the system via 
emission of excitations in the bath (photons, phonons, 
etc) accompanied by transitions between the system en- 
ergy levels. If the coupling is weak, the transition rates 
are small compared to the transferred energy. In the 
classical case, the transitions lead to friction. 

At the quantum level, however, one should take into 
account that the transitions happen at random. The ran- 
domness gives rise to a peculiar quantum noise and the 
related quantum heating of the oscilator. Its important 
manifestation is quantum activation in driven nonlinear 



oscillators [TTl [T2] , where the noise leads to activation- 
type transitions between the states of forced vibrations. 
Quantum activation has been now seen in the experi- 
ment [S]. There is also an observation of the quantum 
heating in the spectrum of a resonantly driven oscillator 
|13j . However, to the best of our knowledge, no direct 
measurements of the relaxation-induced distribution over 
quantum states have been made and no direct means to 
measure this distribution have been proposed. 

In this paper we show that the distribution over the 
states of a modulated nonlinear oscillator can be mea- 
sured spectroscopically. If the oscillator is strongly un- 
derdamped, the power spectrum of its fluctuations and 
the spectrum of the response to an additional weak field 
can display a fine structure. The intensities of the fine- 
structure lines are directly related to the occupation of 
the oscillator states, and the line shapes depend on the ef- 
fective quantum temperature as well as the temperature 
of the thermal reservoir. We note that spectroscopy has 
been long recognized as a means of getting an insight into 
the dynamics of a strongly driven oscillator and, more re- 
cently, of using the oscillator for quantum measurements 
p4l421j . However, the fine structure of the spectra has 
not been discussed earlier. 

We will study the fine structure for an oscillator para- 
metrically modulated at frequency up close to twice the 
eigenfrequency wq. Classically, as a result of parametric 
resonance the modulated oscillator can start vibrating at 
frequency u)p/2. The steady vibrational states are de- 
termined by the balance between the modulation, the 
dissipation due to coupling to a thermal bath, and the 
oscillator nonlinearity [22] • For not too strong modula- 
tion, the oscillator has two vibrational states, which have 
the same amplitude and differ in phase by tt. 

Quantum mechanically, a parametrically modulated 
oscillator can be naturally described in terms of the Flo- 
quet, or quasienergy states ipe{t), such that V'6(i + 2rF) = 



2 



exp{—2ieTp/ti)ipg{t), where Tp = 2Tr/oJF is the modula- 
tion period and e is the quasienergy. The quasienergy 
levels are sketched in Fig. [T] The lowest state in the 
figure corresponds to the stable state of parametrically 
excited vibrations, in the presence of weak coupling to 
the bath; there are two such states in the neglect of tun- 
neling. If these were energy levels, for zero bath temper- 
ature T the oscillator would be at the lowest level and 
its power spectrum would have one line that corresponds 
to the transition to the next level. However, because of 
quantum heating the oscillator occupies higher-£ states in 
Fig. [T] Therefore the spectrum has lines that correspond 
to transitions from these states as well, and the transi- 
tions not only up, but also down in e. Where different 
transitions are spectrally resolved, one gets direct infor- 
mation of the occupation of the states, and thus of the 
quantum temperature that characterizes the quasienergy 
distribution. We note that this temperature generally 
depends on quasienergy [TTl [T^], but for a few lowest 
quasienergy states this dependence is not important. 




FIG. 1. The quasienergy Hamiltonian g near its minimum as 
a function of the coordinate in the rotating frame Q for P — 0. 
The horizontal lines show schematically the quasienergy lev- 
els, which are weakly nonequidistant near the minimum of g. 
The vertical arrows show transitions that lead to the peaks in 
the power spectra with relative intensities that depend on the 
effective Planck number fie. The transitions between higher 
neighboring levels have close, but not identical frequencies, 
which can result in the onset of the fine structure of the noise 
spectrum. In the inset, the arrows indicate the change of g in 
relaxation to and in diffusion away from the classically stable 
state. 



II. OSCILLATOR DYNAMICS IN SLOW TIME 

A. Quasienergy Hamiltonian in the rotating wave 
approximation 

The Hamiltonian of a parametrically modulated non- 
linear oscillator is 

Ho - + [ujI + F cosiujFt)] + Uq\ (1) 



Here, the mass is set equal to one, F is the modula- 
tion amplitude, and 7 is the anharmonicity parameter. 
We assume the modulation to be resonant and compar- 
atively weak, \F\ ^ Wq, so that the nonlinearity is also 
weak, which allows us to keep the lowest-order relevant 
nonlinear term in H, 



\ujf — 2a;o| < 1^0 



l7('z')l«^o'- 



For concreteness we set i^, 7 > 0. 

Following [f 2 , we change to the rotat- 
ing frame using the canonical transformation 
U{t) — exp i^—ia^autFt/'^) , where and a are 
the raising and lowering operators of the oscillator, 
and introduce slowly varying in time dimension- 
less coordinate Q and momentum P, W{t)qU{t) ~ 
Cpar [P cos(wFi/2) " Q sm{oJFt/2)] , C/t {t)pU{t) 



-(CparWF/2) [Psin{ujFt/2) + Qcos{uJFt/2)], 



where 



(2F/37)^/^. The commutation relation between 



P and Q has the form 



[P,Q]^-iX, \^3-fh/FujF. 



(2) 



The dimensionless parameter A plays the role of the 
Planck constant in the quantum dynamics in the rotating 
frame. Respectively, the oscillator raising and lowering 
operators are expressed in terms of P, Q, and A in a stan- 
dard way, 

U^t)aU{t) = (2A)"^/2(P- iQ)exp(-iw_Ft/2). 

It is convenient to analyze the dynamics of a res- 
onantly modulated weakly nonlinear oscillator in the 
rotating wave approximation (RWA). In this approxi- 
mation the Hamiltonian in the rotating frame becomes 
Ho WHoU - ihWU w (FV67) g, where 



g{Q,P) = \{P' 



1 



il-n)P' 



1 



_ ujf{^f - 2tJo) 



The dimensionless operator g describes the oscillator dy- 
namics in slow dimensionless time, with the Schrodinger 
equation of the form of 

= iXdrip — gijj, T = Ft/2ujF- 
The eigenvalues gn of g give the oscillator quasienergies 

En = {F^/&l)gn- 

Operator g does not have the form of a sum of kinetic 
and potential energies and depends on one dimensionless 
parameter ^. We will consider region — 1 < /i < 1 where 
function .g(Q, P) has two minima and a maximum. Its 
cross-section by the plane P = is shown in the inset of 
Fig.[l] 

The full oscillator dynamics, including dissipation, can 
be described by the master equation for the oscillator 
density matrix p. For dissipation that comes from weak 
coupling to a thermal reservoir, which is linear in the os- 
cillator coordinate and possibly momentum, for almost 
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resonant modulation, \2ujo — ujf\ wf, comparatively 
weak nonlinearity, I7I (9^) <C Wq, and the density of states 
of the reservoir weighted with the interaction smooth 
around luq this equation in the RWA has the form 

p= i\^^[p, g] — kp, Hp ^ fi{h + l){a^ ap 
—2apa) + pa} a) + Kn{aa^ p — 2a} pa + paa^). (4) 

Here, operator kp describes dissipation. The transition 
to the interaction representation with respect to the os- 
cillator variables is done using operators U{t), U'^{t). The 
renormalization of the oscillator frequency due to the 
bath is incorporated into wq. The dimensionless param- 
eter K = 2ujpT /F is proportional to the oscillator decay 
rate F; this rate gives the ring-down time 1/2F and the 
quality factor wq/SF, which we assume to be large. We 
note that in Ref. [12] we used rj instead of n. 

In the limit of small n the minima of P) cor- 
respond to the stable stationary states in the rotating 
frame, and thus to the stable states of period-two vibra- 
tions at frequency ujf/2, in the laboratory frame. 



of period-two vibrations. The coordinates of these state 
i(Qoi ^o) in the rotating frame for arbitrary dimension- 
less decay rate k are determined in Sec. |V)seeEq. ([21]); m 
the limit of small k the states are located at the minima 
of g{Q,P), with Qo « (1 + M)^/^fo ~ from Eq. (jsj). 
The states are symmetrical, since they correspond to 
time translation by the modulation period, in the labora- 
tory frame. The contributions to <f>(aj) from fluctuations 
about them are equal, and it is sufficient to study one of 
them. In doing so we will assume that the oscillator is 
localized in the vicinity of the stable state (Qoi^o) ^nd 
disregard interstate transitions. The corresponding term 
in is $0(1^), with 



/•oo 

^o{uj) = Re / dte'- 
Jo 



'{{Sa{t)SaHO))) 



(6) 



Here, 5a{t) — a{t) — ao{t) is the operator a counted off 
from its expectation value oq at the stable state (Qo, Po), 

ao{t) = (2A)-i/2(Po - iQo) exp{-iiJFt/2). 



B. Noise power spectrum 

Of significant interest for the experiment are spectra of 
modulated oscillator [HI [15] , including the power spec- 
trum and the spectra of absorption of an additional weak 
field and radiation emission. Measurements of the power 
spectrum have been already reported f6j for a microwave 
cavity with length effectively modulated by a supercon- 
ducting interference device [l^ ; the related spectrum can 
be studied also through sideband absorption of a Joseph- 
son junction based qubit coupled to a driven nonlinear 
resonator [13]. The power spectrum of the oscillator also 
determines relaxation of a qubit coupled to it [T5] . 

We will consider the power spectrum at frequencies 
close to the eigenfrequency of the parametrically mod- 
ulated oscillator. In the vicinity of the maximum, this 
spectrum is given by 



$H = Re / dte'^'{{a{t)a\0))). 



(5) 



Here, 



{{A{t)Bm) = 



Ijjp 



Air /ujp 



dU{A{t + ti)B{U)), 



where (. . .) indicates ensemble averaging. 

For small fiuctuation intensity, Eq. ( 25 ) , the spectrum 



Q{uj) has distinct peaks near ujf/2. One or two of them, 
possibly with fine structure, see below, come from small- 
amplitude fiuctuations about the stable states. In ad- 
dition, there is an extremely narrow spectral peak from 
rare transitions between the states centered at frequency 



(jJf/2. This peak is discussed in Sec. VI 

Here and in Sees. IllUfVl we will be interested in the 
peak(s) of due to small-amplitude quantum and 

classical fiuctuations about the classically-stable states 



III. QUANTUM TEMPERATURE IN THE 
SMALL DAMPING LIMIT 



A. The Bogoliubov transformation and the 
quasienergy spectrum 

Quantum noise is most clearly manifested in the spec- 
trum if the oscillator relaxation rate is small, so that 
the relaxation-induced width of the quasi-energy levels 
exceeds the inter-level distance. For small A the dimen- 
sionless interlevel distance is \v{g)^ where v{g) is the 
dimensionless frequency of classical vibrations described 
by equations Q = dpg, P — —dqg (in Ref. [T2| we used 
uj{g) instead of i'{g))- The dimensionless level width is 
proportional to the decay rate k. Therefore the condition 
of well-separated levels is v{g) <C k,. 

We are interested in the levels close to the minima of 
function g{Q,P), see Fig. [l] i.e., for g « g^in = ^(1 + 
/i)^/4. Then, from Eq. (|3| the condition of narrow levels 
has the form 

J/O > K, I/O = Hgrain) = 2(1 + p)^/'^ . (7) 

We assume that, at the same time, the level width largely 
exceeds the splitting due to resonant tunneling between 
the minima of g; this splitting is exponentially small for 
A< 1. 

Where these conditions are held, the oscillator motion 
near gmin is weakly damped vibrations at dimensionless 
frequency ~ J^o- It can be studied using the Bogoliubov 
transformation from a, to new operators b, b\ 

U\t)aU{t) = ao{t) + {ub + ^;6^)e~'■"^^*/^ 
. = -^(2.or/^ (l + I) , 

. = -^(2.orV2(i_|). (8) 
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The coefRcients u, v are chosen so that, to second order in 
P,Q — Qo, 9 ~ Xi'ob^b+ const, that is, near its minimum 
g becomes the Hamiltonian of an auxiUary harmonic os- 
cillator with dimensionless frequency i^q. Operators b and 
b^ are, respectively, the lowering and raising operators for 
this auxiliary oscillator. We emphasize that vibrations of 
this oscillator occur in the rotating frame and correspond 
to the vibrations of the original oscillator at frequencies 
ujp/2 ± {F/2ujf)vq- We note that the Bogoliubov trans- 
formation can be written as a squeezing transformation. 



{2\)-^/^{Q-Qn+iP) = fecoshr, - b^ sinhi 



(9) 



with coshr* = iu and sinhr.^ = —iv. 

Higher-order terms in P,Q — Qq in g lead to anhar- 
monicity of vibrations about gmin- In turn, the anhar- 
monicity leads to nonequidistance of the vibrational en- 
ergy levels, that is, of the quasienergy levels of the orig- 
inal oscillator. To the lowest order, the nonequidistance 
is determined by the terms quadratic in b^b taken to the 
first order and by the cubic terms in b, b^ taken to the 
second order. This gives for the eigenvalues of g 



gn « XC'on + ^^^y^^ + 5r 
Here, Vq = vq + XV /2 and 



V 



Ai+ 1 



(10) 



5n 



A. Parame- 
ter V gives the nonequidistance of the levels, with the 
transition frequencies forming a ladder, i^(5„) = (ffn+i — 
gn)IX = vq + XV {n -t- 1) . We note that the frequency step 
XV is proportional to the anharmonicity parameter 7 of 
the original oscillator. Equation ( 10 1 applies for small A 
and small n, where A|T^|n ^ vq. 



B. Master equation in terms of the transformed 
operators 

The full oscillator dynamics near the minima of 



g{Q,P) can be described by the master equation (19) 
with , a written in terms of the operators b^ , b. For 
small K the master equation can be simplified by noting 
that, for K — 0, matrix elements of p on the eigenfunc- 
tions \n) of 5 oscillate in dimensionless time as p„m oc 
exp[— wo(to — n)r/A], for small n,m. Dissipation couples 
matrix elements pmn with Pm'n'- For weak dissipation, 
where Eq. ([t]) holds, the coupling is particularly strong if 
m — n = m—n'. Such coupling is described if in kp writ- 
ten in terms of b\b we keep terms with equal numbers 
of b and W operators, whereas the terms that contain b^ 
and (6^)^ are disregarded. The resulting expression has 
the form 

kp= K{ne + l){b'^bp - 2bpb^ + pb^b) 

{bb''p~2b^pb + pbb^). (11) 

fie = n + (1 + 2n) sinh^ 

= [{p + 2){2n + l)-i^o]/2iyo- (12) 



with 



From the comparison of Eqs. Q (11) one can see that 
fie plays the role of the effective Planck number for vi- 
brations about gmin- 

The stationary solution of the master equation near 
the chosen minimum of g has the form of the Boltzmann 

distribution, pmm oc [fie/{ne + 1)]™ or, in the operator 
form. 



= (fie + ly 

Te = A;/o/ln[(ne 



cxp(-Aj/o6'''6/7^), 
hl)lne]. 



(13) 



Here, Te is the dimensionless effective temperature of vi- 
brations about gmin- Equation (12 1 coincides with the 
results [12] obtained by a completely different method. 
For n = the result coincides also with what follows 
from the analysis of a different model of a modulated os- 
cillator f24| if one uses the appropriate value of the 
squeezing transformation ([9]) . The normalization of p*-***-* 
corresponds to the assumption that the oscillator is lo- 
calized in the vicinity of the stable state (Qo,^o)- The 
distribution over quasienergy states for other systems and 
other relaxation mechanisms were discussed recently in 
Refs.lMlandlMl 



It follows from Eq. ( 12 1 that the effective Planck num- 



ber, and thus also 7^, remain nonzero even for zero tem- 
perature of the bath, fi — T ~ Q. This is a consequence 
of quantum fluctuations that accompany oscillator relax- 
ation, and therefore we call Te quantum temperature. 
The dependence of fie on the dynamical parameter p is 
shown in Fig. [2j It is nonmonotonic, with a minimum at 
exact resonance between the driving frequency and twice 
the oscillator eigenfrequency, where p oc up — 2wo — 0. 
The value of fie increases with decreasing p + 1 oc i^q, 
i.e., close to the bifurcation point where the period- two 
vibrations are excited. However, the assumption j/q S> k 
breaks down sufficiently close to the bifurcation point, 
which imposes a restriction on fie- 

The occurrence of the minimum of fie is an interesting 
feature of the parametrically modulated oscillator. For 
p — 0, from Eq. (12) fie — fi] the effective temperature 



is equal to the thermal bath temperature, and 7^ = for 
T — Q. We note that this is an asymptotic result that 
applies only very close to (?niin- The stationary distribu- 
tion pirn is Boltzmann-like only to the leading order in 
the distance gn — gmin <C (7min- Strictly speaking, the ef- 
fective temperature is quasienergy-dependent [12] . Even 
for p = T = Q quasienergy states with n > 1 are oc- 
cupied, but in the range of small n this occupation is 
much smaller ioi p — {) than for \p\ ^ 1. The role of the 
corrections to the Boltzmann distribution for /i = is 
illustrated in the inset of Fig. [2] 
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FIG. 2. The efTective oscillator Planck number Tie SI'S ^ func- 
tion of the scaled frequency detuning fi — up {up — 2u)o)/ F for 
different values of the Planck number n, for small decay rate. 
The inset shows the populations of the low-lying quasienergy 
states pnn for = 0. The solid lines in the inset show the 
results of the calculation that takes into account relaxation- 
induced transitions between neighboring quasienergy levels, 
Eq. (121, shown in the main figure. The dots, squares, and 



circles show the result of the full calculation, Ref. J^. 



IV. FINE STRUCTURE OF THE POWER 
SPECTRUM 

A. General expression for the spectrum near its 
maximum 



Equations (10 1 - (11) describe the dynamics of the 
modulated underdamped oscillator in terms of an aux- 
iliary oscillator in thermal equilibrium with temperature 
Te- The power spectrum <i>o(<^) of the original oscillator 
near its maxima can be expressed in terms of the power 
spectrum of the auxiliary oscillator. From Eqs. (|6| and 
([S]), for frequencies w — ujp/2 close to the frequency of 
vibrations about the stable states. 



20; J 



\u\^^b{v), (14) 



/•OO 

<i>b(i^) = Rc / dTe--((6(T)&t(0)))^„t. 
Jo 

Here, v is the dimensionless frequency counted off from 
lof/2, 

V = {2ujf/F) {w - w_f/2) ; 

the subscript in ((. . .))rot indicates that the correlator is 
calculated in the rotating frame, 

((A(r)S(0))),ot = Tr Ap(r; B), p(0; B) = Bp<^''\ 

(15) 

where p(t; B) satisfies the master equation Q with the 
dissipative term of the form Eq. (11) and p^^*) is the 
stationary distribution given by Eq. (13). 



In deriving Eq. (14) we took into account that, if 



about 5min, exp(i gr/A)^ exp(— i gr/ A) w exp(— wo''')^, 
and therefore function $&(z^) describes the dominating 
contribution to $o for close to i^q. Using the fact that 
small-amplitude vibrations near gmin can be thought of 
as being close to equilibrium, one can show that the 
peak of $0 for i' close to —i/q is described by function 
{2ujp\v\yF) exp(-Ai.o/re)$;,(-i^). 



B. Effective partial spectra representation 

The problem of the power spectrum of a weakly non- 
linear underdamped oscillator was discussed previously 
[27] . Applying the results to the spectrum $b(i^) of the 
auxiliary oscillator, after some straightforward transfor- 
mations we obtain 



2$b(:/) (fie + l)Re Mn, v)\ 

n 

(j)b{n, v) = 4n(A - l)"-i(A + l)-("+i) 
X [K(2Hn -I -id)- i{v - V(i)y^ 



(16) 



where 



A = N^^l + ii9(2ne + 1)] , d 



= \V/2k, 
(Re H > 0). 



Equation ( 16 ) can be thought of as a representation 
of the spectrum as a sum of effective partial spectra 
Re 4>b{n,v) that correspond to transitions n — 1 — >■ n 
between the oscillator quasienergy levels. Functions 
(/)„(n, I/) depend on two parameters, d and n^.. Param- 
eter gives the ratio of the difference XV = v{gn) — 
^{gn-i) between neighboring transition frequencies and 
the broadening k of the quasienergy levels, whereas the 
effective Planck number fie gives the typical width of the 
stationary distribution over the levels. 

The form of (f>b{n-, v) is particularly simple for a com- 
paratively large frequency spacing or small damping, 
A|y| ^ K,. Note that for small A this is a much stronger 
restriction on the decay rate than the condition n <^ 



used to derive Eq. (16). For such a small decay rate 
(t)b{n,u) « _ " ( 

Tip + 1 



^-\uo[n-l)/T^ 



we ignore dissipation and the nonlinearity of vibrations 



{Kr,-i[v-v{gn-l)]} \ > 1, 

K„ = K[2n(2ne + 1) - 1]. (17) 

In this limit Re (j)b{n,h') is a Lorentzian line centered at 
the frequency ly^gn-i) = vo+XVn of transition n—l — >■ n, 
with halfwidth k„ equal to the half-sum of the reciprocal 
lifetimes of the levels n—l and n. Associating (j)b{n,i') 
with a partial spectrum is fully justified in this limit. 
The amplitude of ^^(n, v) contains the Boltzmann factor 
proportional to the population of the quasienergy level 
n—l. The overall spectrum <I>f,(j^) has a pronounced fine 
structure for \'d\ 3> 1. The line intensities immediately 
give the effective quantum temperature 7^. 
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As |t?| decreases the partial spectra start to overlap and 
for I'd] < 1 they can no longer be identified. Indeed, the 
typical dimensionless time < the oscillator spends 
in a given quasienergy state becomes smaller than the 
distance |AV^|^^ between different transition frequencies 
^{gn). Therefore such frequencies cannot be resolved. 
In the limit |t?| — we have j/) oc (5n,i, and the 

spectrum has the form of a single Lorentzian peak of 
dimensionless width k. 



$t(i/) = (ue + 1) Re [k - i{v - uo)Y 



(18) 




15 / ~ w 



FIG. 3. The fine structure of the noise spectrum of a mod- 
ulated oscillator near resonant frequency uip/i + Fvo/2u}f 
[u/k, = {u) — ujf/'2)/T\. The intensities of the fine structure 
lines are determined by the effective Planck number fie- Due 
to quantum heating, the fine structure is seen even for zero 
bath temperature. The widths of the individual lines increase 
with increasing oscillator Planck number n, leading to the 
smearing of the fine structure. 

The evolution of the spectrum with varying fie depends 
on \{}\. For |t?| ^ 1 the spectrum is Lorentzian for fie ^ 1, 
but becomes non-Lorentzian and asymmetric for large fie, 
where l^^lfie > 1- However, the spectrum does not have a 
fine structure in this case. On the other hand, for ^ 1 
a fine structure emerges, but only in a limited range of n^. 
It is seen from Eq. (17) that for fie ^ 1 only ^{,(1, z/) has 



an appreciable intensity, </'b(n 7^ l,i^) are small. On the 
other hand, for large fie the linewidth k„ becomes large 
and the spectral lines with different n overlap, starting 
with large n. The evolution of the fine structure with 
varying n is illustrated in Fig[3j The quantitative results 
confirm the above qualitative arguments. 



V. MODERATE DAMPING 
A. Master equation in the Wigner representation 

The oscillator power spectrum should be analyzed dif- 
ferently if the oscillator damping is not that small. We 
assume that the original oscillator remains underdamped. 



Fk/u)p <^ ujQ, but for the auxiliary oscillator which vi- 
brates about the minimum of g{Q,P), the dimension- 
less decay rate k exceeds the nonequidistance of the 
quasienergy levels, k 3> A|y|, so that the power spectrum 
does not have the fine structure discussed in Sec. IIIII 
The latter condition indicates that the quantum effects 
related to the difference of the transitions frequencies are 
small, since X (x h. However, other quantum effects are 
still important, as seen below, and in particular the spec- 
trum strongly depends on the quantum temperature. We 
note that, for k 3> A|y|, the ratio of the width of the 
quasienergy levels to the distance between them oc k/vq 
can still be arbitrary. 

The analysis of the power spectrum for moderate 
damping can be done by writing master equation Q in 
the Wigner representation, 

pw = -V (Kpw) + XL^^^pw + X^W'^pw (19) 

Here, pw is the density matrix in the Wigner represen- 
tation, 



Pw{Q.P:t)^ J die-'^P'^p {q 



where p{Qi,Q2',t) = {Qi\p{t)\Q2) is the density matrix 
in the coordinate representation. In Eq ( 19 ) we use vector 



notations, K = {Kq,Kp) and V = {dq.dp). 

Vector K determines the evolution of the density ma- 
trix in the absence of quantum and classical fluctuations, 

Kq = dpg - Kp = -dgg - fiP, (20) 



whereas the terms cx A in Eq . (|1 9[) account for fluctua- 
tions. If we set A 
motion Q — K, 



0, Eq. ( 19 1 will describe classical 
P — Kp. The condition K = gives 



the values of Q, P at the stationary states of the oscillator 
in the rotating frame. For \p\ < (1 — k^)^/^ the system 
has 3 stationary states. One is located at Q = P = and 
is unstable. The other two are located symmetrically at 
±((3o,Po) with Qq — rocos6', Pq = rpsin^, where 



2N1/2 



(21) 



and 9 = arctan { [l — (1 — k^)^/'^] /k}. These states are 
asymptotically stable. Respectively, the real parts of the 
eigenvalues of matrix fC, 



/C„- = [dK,/dXj 



:>:Po 



(Xi = Q, X2 = P) (22) 



are negative [the subscript {Q(),Po) indicates that K. is 
calculated at point {Qq, Pq)]. For k ^ 1 the stable states 
correspond to the minima of g{Q, P) in Fig. [ij 

The terms L^^^ and L'^' in Eq. (19 1 describe, respec- 



tively, the fluctuations due to decay processes and purely 
quantum fluctuations that are not related to the coupling 
to a thermal bath. 



=K(f^+l/2)V^ 
L^^) = -^^{Qdp - PdQ)V^ 



(23) 
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The decay-related fluctuations lead to diffusion in (Q, P)- 
space, as seen from the structure of L^^\ In contrast, the 
term L^^^ is independent of k and contains third deriva- 
tives; for small A it is not important close to the stable 
states. 

It follows from Eqs. (19) and (23) that, in the 



Wigner representation, the stationary distribution p'^^ 
has Gaussian peaks at the stable states ±{Qo,Po). They 
are of the same form for the both states, and close to 
(Qo,Po) 



(deti) 



1/2 



2n + l 
2kA^ + At + t"^ A = 0. 



exp 



(5Xyl(5X 

'A(2n+ 1) 



(24) 



where 



5y.= {5Q,5P) = {Q-Qo,P-Pa) 



is the distance from the stable state, |(5Xp <C Qq 



Pi 



Equation ( 24 ) shows that the condition for quantum and 



classical fluctuations to be small is 



A(2n+ 1) < rl. 



(25) 



From Eqs. ( 20 ) , ( 22 ) , and ( 24 ) Tr A = 2 is independent of 



the parameters of the system. However, the distribution 
( 24 ) is squeezed [2H] , and the variances oi Q — Qq and 
P — Po depend on the oscillator parameters. 

It is useful to note that, in the small-damping limit 
K <C i^o, matrix A is diagonal, with All ~ 2(1-1- /i)/(2-|-/i), 
A22 ~ 2/(2 + /i), and A12 oc h/vq. Using the rela- 
tion 2ne -f- 1 = (m + 2)(2fi -f that follows from 



Eq. (12), one can see that the above expression for A 

leads to p'^^ cx exp{-2[9((5, P) - .gmi„]/[Ai/o(2ne + 1)]}, 
which is the standard form of the Wigner distribution 
of a harmonic oscillator; in the present case, the result 



refers to the auxiliary oscillator discussed in Sec. Ill with 
Hamiltonian g{Q,P), frequency uq, and Planck number 
fie- The result is fully consistent with what was found in 



Sec. Ill using a different method. 



B. Power spectrum for moderate damping 



Equations (19) and (24) allow one to flnd the oscil- 



lator power spectrum for an arbitrary relation between 
the width of the quasienergy levels and the level spacing 
h/uq but for X\V\ <C k. Again, we will be interested in 
the contribution to the spectrum from small-amplitude 
fluctuations about the stable states. The general expres- 
sion for this contribution follows from Eqs. ([6]) and (15), 



F 
2lup 



$0(0.)= Re P dre^'^^ j ^^('5P 



i6Q) 



(26) 



x^1^^(Q,P;t), 
where function p\^' satisfies master equation (|19|) with 
the initial condition 



p[^)(Q,P;0) = 2 



5P + i5Q - -\{idQ + dp) 



(27) 



Function p''y/\Q, P; 0) is the Wigner-transform of the op- 
erator ((5P -I- iSQ)p{T — 0); for operator p this trans- 
form is defined by Eq. ([20|. Factor 2 in Eq. ([27]) ac- 
counts for the contribution of fluctuations about the state 
— (Qoi-Fb); we have also taken into account in Eq. (26) 
that J dQdPpw 



27rA. 



The calculation of the power spectrum using Eqs. (26) 
and 



\21\ is similar to that performed in the classical 
|14j and quantum theory |18j for the power spectrum 
of an oscillator modulated by an additive force at fre- 



quency close to loq. One should replace pw in Eq. (19) 
with p'^y^^ , set K sa /S(5X, multiply the equation by 
exp(ii'T) and then in turns by 5P and 5Q. One should 
then integrate the resulting equation over t,P,Q, as in 
Eq. ([26|. This will lead to two coupled linear equations 
for the Fourier transforms of ((5P(t)[(5P(0) + i5Q{0)]) 
and {6Q{t)[SP{0) +iSQ{0)]). The inhomogeneous parts 
of these equations are determined by the average values 
(SXiSXj), which can be found from Eq. (24). A straight- 



forward but cumbersome calculation gives 



F 



2rl 



n(l 



-aV2) 



UJp{2uj — LJp) 

F ' 



(28) 



where tq is the dimensionless amplitude of parametrically 
exci ted vibrations in the neglect of fluctuations given by 
Eq. (21). The frequency Va = 2ro(rQ— /i)^/^ characterizes 



damped vibrations about the stable state in the absence 
of fluctuations, — ~ det JC> Q. 

For small but not too small damping, \\V\ ^ n <^ vq 
we have Va~i^a- One can then show from Eq. (28) that 



function ^q{uj) has two Lorentzian peaks at dimension- 



less frequencies ij/Q with halfwidth k. The expression for 



the peak at t'o coincides with Eqs. ( 14 ) and ( 18 ), whereas 
for V close to —vq 



2wp 



$o(w) 



\V\ UeK 



2 ' ^21-1^ 



(29) 



in agreement with Sec. |IV[ We emphasize that, in the lab- 
oratory frame, the spectral peaks described by Eqs. ( 18 ) 
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and (291 lie on the opposite sides of frequency ujf/'2, at 
the distance Fvq/2ujp in dimensional frequency. The 
ratio of their intensities is proportional to the factor 
exp(— Ai'o/T^) and thus strongly depends on the quantum 
temperature, which provides an independent means for 
measuring this temperature. Even for zero temperature 
of the thermal reservoir 7^ > 0, generally, and therefore 
both peaks are present in the spectrum. 

The evolution of the spectrum ([28]) with varying oscil- 
lator parameters is illustrated in Fig. |4] With increasing 
bath temperature the intensity of the spectrum goes up, 
but the peaks at vq and —vq remain different. The ratio 
of their intensities approaches for n ^ l.On the 

other hand, with decreasing Va/ n the peaks start overlap- 
ping, and ultimately form a single peak. For small Va/i^ 
the peak is centered ai v = Q (at uj = iopjl^ in the lab- 
oratory frame) and has halfwidth v\I1k <C k. The limit 
Val ^ 1 is relevant for the vicinity of the bifurcation 



point /i = — (1 
disappear. 



where the period- two vibrations 




VI. SUPERNARROW SPECTRAL PEAK AND 
SOME GENERALIZATIONS 



o 
CM J 







-1 


1 




V 


(b) n = 


f\K = 0.2 






K = 0.35 / 




Av = 0.5^/^ 





Along with small-amplitude fluctuations around the 
stable vibrational states, quantum and classical fluctu- 
ations lead to occasional interstate switching. Unless the 
damping rate is extraordinarily small, even for zero bath 
temperature the switching occurs via transitions over the 
quasienergy barrier that separates the minima oig{Q, P), 
see Fig.[l]jT^. For small A(2n-|-1) the transition rate Wtr 
is exponentially small, Wtr oc rexp(— i?/A) (the dimen- 
sional decay rate F ~ kF/2ujf)', the effective activation 
energy R was discussed earlier [11] [T^] . 

An important manifestation of interstate switching is 
the occurrence of an additional peak in the oscillator 
power spectrum. It is centered at frequency ujp/2 and is 
supernarrow in the sense that its width is much smaller 
than F. The peak is analogous to the supernarrow peak 
in the spectra of oscillators with coexisting vibrational 
states in a resonant additive field [HI [TB] . The distinc- 
tion is that, for a parametrically modulated oscillator, 
average populations of the stable states are equal for all 
parameter values in the range of bistability, as a conse- 
quence of the symmetry with respect to time translation 
by 2-K /lof- The supernarrow peak in the response of the 
parametric oscillator was discussed earlier , and in the 
noise spectrum of a parametric oscillator it was recently 
seen in the experiment [H]. 

To describe the peak in the noise spectrum we note 
that the populations and p- of the stable vibrational 
states (Qo:^o) and — ((5o,Po), respectively, satisfy the 
balance equation 



dp±/dt = ±Wtr{p-- p+). 



(30) 



FIG. 4. The scaled contribution to the power spectrum 
from small-amplitude fluctuations about the stable period- 
two states [F /2ujf)^o{lli) for moderate damping where the 
fine structure is smeared out. The data refer to /i = —0.75. 
Panels (a) and (b) show the evolution of the spectra with vary- 
ing bath temperature, which determines the Planck number 
n, and with the dimensionless decay rate respectively. For 
comparatively small n the spectrum has two peaks, which 
in the laboratory frame are located at ~ top 12 ± Fva/^uip. 
They correspond to transitions up and down in quasienergy, 
cf. Fig.[l] The peak at liJf /2 — Fva/2u)F for zero bath temper- 
ature emerges due to the quantum heating. With increasing 
decay rates the peaks merge. Superimposed on the shown 
spectra is the supernarrow peak at frequency top 12 (i^ = in 
the rotating frame), which is described by Eq. (33 1. 



a(t),a'f(t) between the stable-states values (ao(t), aQ(t)) 
and — (ao(t), aQ(t)); we note that fluctuations about the 
stable states are averaged out on time scale ~ F~^ <C 
Wt7^. The contribution of these fluctuations to the time 
correlation function of a, is 

((a(t)at (0)))tr « aoe-^""*/'[p+(t; a^) - P-{t; «^)], (31) 



Fluctuations of the populations p± lead to fluctu- 
ations of the expectation values of the operators 



where p±{t;a^) satisfy Eq. (30 1 with initial conditions 
p±{0;a}) = ±ao/2 that follow from the stationary state 
populations being equal to 1/2. 

From Eqs. and ^ we obtain the full ex- 

pression for the power spectrum as 

$(c^) = $o(w) + $tr(c^), (32) 
where $tr describes the interstate-transition induced con- 
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tribution, 



$tr(w) = 2W^tr|aol 



(33) 



Function $tr('^) has the shape of a Lorentzian peak with 
halfwidth 2Wtr <C F. The intensity of this supernar- 
row peak is determined by the squared scaled amphtude 
of the period- two vibrations oc Pq + Qq. The area of 
the peak is independent of the bath temperature, but its 
width sharply increases with the increasing temperature. 



A. Quantum temperature for zero-amplitude states 

In the parameter range > (1 — k^)^/^ the oscilla- 
tor has a stable state where the amplitude of vibrations 
at frequency ujp/2 is zero. Even though the oscillator 
does not vibrate on average, fluctuations about the zero- 
amplitude state are modified by the periodic modulation. 
These fluctuations are described by Eqs. ([3| and Q. For 
small damping, in the frame rotating at frequency wf/2, 
fluctuations are random vibrations of an auxiliary oscilla- 
tor at a dimensionless frequency v'q — (/i^ — 1)^/^. These 
vibrations can be described using the Bogoliubov trans- 
formation similar to that in Sec. |III[ with — Q and 
with u and v replaced by u' and v' , respectively, 



|l/2 



|l/2 



'1/2 



(I 



1 



1/2 



II 



(34) 



The effective Planck number of the vibrations of the 
auxiliary oscillator is 



2n + l 



II 



Ml 



|i-mI)- 



1 



(35) 



Even where n = 0, we have a nonzero n'^ — |u'p > 
0. The effective temperature of the auxiliary oscillator 
increases close to the bifurcation points /i « ±1 where 
the zero-amplitude states of the original oscillator loose 
stability, with n'^ w |/i2-l|-i/2/2 for fi = 0. On the other 
hand, far from the bifurcation points, where /i^ 3> 1, we 
have n'g ss n, i.e., as expected for a zero-amplitude state, 
the temperature of the auxiliary oscillator approaches the 
bath temperature. 

The distribution over the quasienergy states for the 
oscillator fluctuating about a zero-amplitude state can 
be directly measured spectroscopically through the fine 
structure of the power spectrum. The nonequidistance 
of the quasienergy levels in this case is the same as for 
the energy levels of the original oscillator in the absence 
of modulation. 



VII. CONCLUSIONS 

This paper was focused on quantum fluctuations that 
accompany relaxation in modulated oscillators. These 



fluctuations lead to a finite width of the distribution of 
the oscillator over quasienergy states, even for zero tem- 
perature of the thermal bath that causes relaxation. We 
call this effect quantum heating. It gives an extra con- 
tribution to the standard quantum fluctuations related 
to a finite width of the oscillator distribution over the 
coordinate and momentum in each quasienergy state. 

As a consequence of the finite width of the quasienergy 
distribution, the power spectrum of an underdamped os- 
cillator in the rotating frame has peaks at frequencies 
that correspond to transitions with increasing or decreas- 
ing quasienergy. These frequencies have opposite signs, 
and the ratio of the peak amplitudes is determined by 
the width of the quasienergy distribution. The transi- 
tions occur primarily between neighboring quasienergy 
levels. We note that, in the language of quantum optics, 
one can think of the peaks as resulting from parametric 
down-conversion: a photon at frequency up splits into 
photons at frequencies up /2±5uj. However, the processes 
involved are substantially multiphoton, and a description 
in terms of quasienergies is more adequate. 

We have shown that the peaks of the power spectra 
may have fine structure. It emerges where the difference 
in frequencies of transitions between neighboring pairs 
of quasienergy levels exceeds the decay rate. In dimen- 
sionless units this condition has the form X\V\ ^ k, see 



Sec. IV The power spectrum of a nonlinear oscillator 



may display fine structure also in the absence of periodic 
modulation, provided the nonequidistance of the energy 
levels exceeds their width, which is the same condition 
but applied to the energy rather than quasienergy levels. 
Quantitatively, it has the form A ^ 2/c [27] . 

For the period- two states \V\ > 2, and \V\ becomes 
large near the bifurcation point where /i + 1 is small, see 
Eq. (10). Hence it can be significantly easier to observe 



the fine structure for a modulated oscillator than for an 
unmodulated one. In addition, the observation does not 
require that the excited states of the unmodulated os- 
cillator be thermally populated. A comparatively strong 
nonequidistance of the energy levels of unmodulated os- 
cillators has been already achieved in circuit QED; in par- 
ticular, it underlies the operation of the transmon qubits 
|30 j . Therefore the fine structure predicted in this paper 
should be accessible to the experiment. A similar fine 
structure can be observed in the power spectrum of an 
oscillator driven by an additive force with frequency close 
to the oscillator eigenfrequency. 

An interesting feature of the parametrically modulated 
oscillator, that does not occur in an additively driven os- 
cillator, is the occurrence of the parameter value where 
the effective temperature of the quasienergy distribution 
coincides with the temperature of the bath, to the lead- 
ing order in the distance from the stable state along the 
quasienergy axis. Near bifurcation points where the sta- 
ble state disappears, on the other hand, the effective tem- 
perature sharply increases. 

Another important feature is the supernarrow peak at 
frequency up/2, which emerges in the response [29) and 
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also in the noise spectrum, where it has been aheady seen 
in the experiment [6]. In contrast to the supernarrow 
peak for additively driven oscillators [TH [TB] , for para- 
metric oscillators the peak has large intensity in a broad 
parameter range, everywhere where the period-two states 
are significantly populated. The width of the peak is de- 
termined by the rate of switching between the period-two 
states and is much smaller than the oscillator relaxation 
rate. 

For small damping, k, <^ the quantum- heating 

induced fine structure should be observable not only in 
the noise spectrum, but also in the spectrum of lin- 
ear response to an additional weak field at frequency uj 
close to ujp /2 ± Fvq/2ujp. Since near its maximum the 
quasienergy distribution is of the Boltzmann form, this 
spectrum can be analyzed using an appropriately modi- 
fied fluctuation-dissipation relation. Its shape is similar 
to that described by Eq. ( 16 ). Where the fine structure is 



smeared out, k ^ -^l^li quantum effects weakly change 
the response to an extra weak field, and the analysis of 
this response for a parametrically modulated oscillator 
can be done in the same way as for a classical oscilla- 
tor driven by a resonant additive force [H]. We note 



that the response emerges both at frequency uj and the 
mirror frequency uip — uJ- 

In conclusion, we have demonstrated that quantum 
fluctuations, which accompany relaxation of a period- 
ically modulated oscillator, can be observed by study- 
ing the oscillator power spectrum. For a parametrically 
modulated oscillator, we found the spectrum in an ex- 
plicit form. In the laboratory frame, the spectrum may 
have two peaks located on the opposite sides of half the 
modulation frequency, or, for higher damping, a single 
peak. Where the spectrum has two peaks, the ratio 
of their intensities is determined by the quantum tem- 
perature, which characterizes the distribution over the 
quasi-energy states of a modulated system. Generally, it 
exceeds the bath temperature. For small damping, the 
spectral peaks may display a fine structure. The inten- 
sities of the fine-structure lines as well as their shapes 
are also determined by, and sensitively depend on the 
quantum temperature, suggesting an independent way 
of measuring it. 
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